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ABSTRACT:  For  nonlinear  functions  f  of  a  random  vector  Y,  E[f(Y)]  and  f(E[Y])  usually  differ.  Conse¬ 
quently  the  mathematical  expectation  of  Y  is  not  intrinsic:  when  we  change  coordinate  systems,  it  is 
not  invariant. This  article  is  about  a  fundamental  and  hitherto  neglected  property  of  random  vectors  of 
the  form  Y  =  f  (X(f) )  ,  where  X(f)  is  the  value  at  time  t  of  a  diffusion  process  X:  namely  that  there 
exists  a  measure  of  location,  called  the  "intrinsic  location  parameter"  (ILP),  which  coincides  with 
mathematical  expectation  only  in  special  cases,  and  which  is  invariant  under  change  of  coordinate 
systems.  The  construction  uses  martingales  with  respect  to  the  intrinsic  geometry  of  diffusion  pro¬ 
cesses,  and  the  heat  flow  of  harmonic  mappings.  We  compute  formulas  which  could  be  useful  to  stat¬ 
isticians,  engineers,  and  others  who  use  diffusion  process  models;  these  have  immediate  application, 
discussed  in  a  separate  article,  to  the  construction  of  an  intrinsic  nonlinear  analog  to  the  Kalman  Fil¬ 
ter.  We  present  here  a  numerical  simulation  of  a  nonlinear  SDE,  showing  how  well  the  ILP  formula 
tracks  the  mean  of  the  SDE  for  a  Euclidean  geometry. 

RESUME:  Pour  une  fonction  non  lineaire  f  d'un  vecteur  aleatoire,  E[f(Y)]  et  f(E[Y])  sont 
usuellement  differents.  Par  consequent,  I'esperance  mathematique  de  Y  n'est  pas  intrinseque:  quand 
nous  changeons  le  systeme  des  coordonnees,  elle  n'est  pas  invariante.  Cet  article  concerne  une  pro¬ 
priety  fondamentale,  negligee  jusqu'a  maintenant,  des  vecteurs  aleatoires  de  la  forme  Y  =  f  (X  (t) )  , 
ou  X(f)  est  la  valeur  au  temps  t  d'un  processus  de  diffusion  X:  c'est  a  dire  qu'il  existe  une  mesure  de 
position,  nommee  le  "parametre  intrinseque  de  centrage"  (PIC),  qui  coincide  avec  I'esperance 
mathematique  seulement  dans  des  cas  specifiques,  et  qui  est  invariante  par  changement  du  systeme 
des  coordonnees.  La  construction  utilise  des  martingales  en  rapport  avec  la  geometrie  intrinseque  des 
processus  de  diffusion,  et  le  flot  de  chaleur  des  applications  harmoniques.  Nous  calculons  des  for- 
mules  qui  peuvent  etre  utiles  aux  statisticiens,  aux  ingenieurs,  et  a  toute  autre  personne  qui  utilise  des 
modeles  fondes  sur  des  processus  de  diffusion;  ces  formules  se  mettent  en  service  a  la  construction 
d'une  analogue  non  lineaire  intrinseque  du  filtre  de  Kalman,  discutee  dans  un  autre  article.  Nous 
presentons  ici  une  simulation  numerique  d'une  EDS  non  lineaire,  qui  montre  la  precision  avec  laque- 
lle  la  formula  de  PIC  suit  la  moyenne  de  I'EDS  pour  une  geometrie  Euclideenne. 
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Technical  Overview 


1  Technical  Overview 

Suppose  X  is  a  Markov  diffusion  process  on  ,  or  more  generally  on  a  manifold  N.  The  diffusion 
variance  of  X  induces  a  semi-definite  metric  (.|.)  on  the  cotangent  bundle,  a  version  of  the  Levi-Civita 
connection  F,  and  a  Laplace-Beltrami  operator  A.  We  may  treat  X  as  a  diffusion  on  N  with  generator 
(1/2)  A,  where  ^  is  a  vector  field. 

For  sufficiently  small  5  >  0,  Xg  has  an  "intrinsic  location  parameter",  defined  to  be  the  non-random 
initial  value  Vq  of  a  F-martingale  V  terminating  at  Xg .  It  is  obtained  by  solving  a  system  of  forward- 
backwards  stochastic  differential  equations  (FBSDE):  a  forward  equation  forX,  and  a  backwards 
equation  for  V.  This  FBSDE  is  the  stochastic  equivalent  of  the  heat  equation  (with  drift  for  harmonic 
mappings,  a  well-known  system  of  quasilinear  PDE. 

Let  N,t>0}  be  the  flow  of  the  vector  field  and  let  =  (xq)  e  N  .  Our  main  result  is 

that  exp  V„  can  be  intrinsically  approximated  to  first  order  in  T  N  by 

Xg  U  Xg 

Vd(^g(xo)  (Hg)  -j^((^g_p,(Vd(^^(xo))dn^ 

r* 

where  n  =  ((|)  ).(.|.)  ds  e  T  N®T  N  .  This  is  computed  in  local  coordinates.  More  generally, 

^  ^  ^0  Xg  2 

we  find  an  intrinsic  location  parameter  for  \|/  (Xg)  ,  if  i|/:N  M  is  a  C  map  into  a  Riemannian 
manifold  M.  We  also  treat  the  case  where  Xq  is  random. 

2  Geometry  Induced  by  a  Diffusion  Process 


2.1  Diffusion  Process  Model 

Consider  a  Markov  diffusion  process  {X^,  t  >  0}  with  values  in  a  connected  manifold  N  of  dimension 
p,  represented  in  coordinates  by 


dx[  =  b'{X^)dt+  Y^a]{X^)dW^, 


where  is  a  vector  field  on  N,  a  (x)  =  {o'-  (x))  e  L  ,  and  W  is  a  Wiener  process  in 

D  c/x  •  i  i  '2. 

R^  .We  assunie  for  simplicity  that  the  coefficients  b  ,  a.  are  C  with  bounded  first  derivative. 


2.2  The  diffusion  variance  semi-definite  metric 

Given  a  stochastic  differential  equation  of  the  form  (1)  in  each  chart,  it  is  well  known  that  one  may 

2 

define  a  C  semi-definite  metric  (.|.)  on  the  cotangent  bundle,  which  we  call  the  diffusion  variance 
semi-definite  metric,  by  the  formula 


(2) 
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Note  that  (.|.)  may  be  degenerate.  This  semi-definite  metric  is  actually  intrinsic:  changing  coordinates 
for  the  diffusion  will  give  a  different  matrix  (a.)  ,  but  the  same  semi-definite  metric.  We  postulate: 

Axiom  A:  The  appropriate  metric  for  the  study  of  X  is  the  diffusion  variance  semi-definite  metric,  not 
the  Euclidean  metric. 

The  p  X  p  matrix  ( (a  •  a)  defined  above  induces  a  linear  transformation  a  (x)  :T^*N  T^N , 
namely 

a  (x)  (dx)  =Y^{a  ■  a)  ''d/dx. . 

Let  us  make  a  constant-rank  assumption,  i.e.  that  there  exists  a  rank  r  vector  bundle  E  N ,  a  sub¬ 
bundle  of  the  tangent  bundle,  such  that  =  range  (a  (x))  cT^N  for  all  xe  N  .  In  Section  7  below, 
we  present  a  global  geometric  construction  of  what  we  call  the  canonical  sub-Riemannian  connec¬ 
tion  V°  for  (.|.),  with  respect  to  a  generalized  inverse  g,  i.e.  a  vector  bundle  isomorphism 
g:TN^T*N  such  that 

a(x)  •gix)  •afx)  =  a(x)  .  (3) 

In  local  coordinates,  g  (x)  is  expressed  by  a  Riemannian  metric  tensor  (g,.^)  ,  such  that  if 
=  (a  ■  a)  \  then 


The  Christoffel  symbols  {rf}  for  the  canonical  sub-Riemannian  connection  are  specified  by  (83) 
below.  The  corresponding  local  connector  F  (x)  e  L  (T^R^  ®  fT^R^)  can  be  written  in  the  more 
compact  notation: 


2g(r(x)  (u®v))  •  w  =  D<g  (v)|g  (w))(u)  H- D<g  (w)  |g  (u))  (v)  -  D<g  (u)  |g  (v))  (w)  ,  (5) 

where  g  (F  (x)  (u  ®  v) )  is  a  1  -form,  acting  on  the  tangent  vector  w. 


2.3  Intrinsic  Description  of  the  Process 

The  intrinsic  version  of  (1 )  is  to  describe  X  as  a  diffusion  process  on  the  manifold  N  with  generator 

L  =  ^  +  (6) 

where  A  is  the  (possibly  degenerate)  Laplace-Beltrami  operator  associated  with  the  diffusion  variance, 
and  ^  is  a  vector  field,  whose  expressions  in  the  local  coordinate  system  {x\  ...,x^}  are  as  follows: 


(7) 
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G-Martingales 


3  r-Martingales 

-  2 

Let  r  be  a  connection  on  a  manifold  M.  An  H  F-martingale  is  a  kind  of  continuous  semimartingale 
on  M  which  generalizes  the  notion  of  continuous  martingale  on  :  see  Emery  [1 5]  and  Darling 
[7].  We  summarize  the  main  ideas,  using  global  coordinates  for  simplicity. 

Among  continuous  semimartingales  in  Ito's  formula  shows  that  local  martingales  are  character¬ 
ized  by 

f(X^)  -f(Xo)  -  (l/2)f^D^f(X^)  {dX®dX)^e  Vfe  C^R")  ,  (8) 

where  (dX  ®  dX) ''  is  the  differential  of  the  joint  quadratic  variation  process  of  X'  and  X^ ,  and 
refers  to  the  space  of  real-valued,  continuous  local  martingales  (see  Revuz  and  Yor  [24]).  For  vector 
fields  ^  on  R'^,  and  co  e  (R'^)  ,  the  smooth  one-forms,  a  connection  F  gives  an  intrinsic  way  of 
differentiating  co  along  ^  to  obtain 

V^coe  Li^R"')  . 

V^co  •  ^  is  also  written  V(o(^®Q  .When  co  =  df,  this  gives  the  Hessian 

Vdf(Dj@D.)  = 

/ 


E  <X,  X)^  =  E]^ ^g,.  (X^)  d<X'  X)t  <  -  .  (11) 

The  F-martingale  Dirichlet  problem,  which  has  been  studied  by,  among  others,  Emery  [14],  Kendall 
[18],  [19],  Picard  [22],  [23],  Arnaudon  [1],  Darling  [6],  [7],  and  Thalmaier  [27],  [28],  is  to  construct 
a  F-martingale,  adapted  to  a  given  filtration,  and  with  a  given  terminal  value;  for  the  Euclidean  con¬ 
nection  this  is  achieved  simply  by  taking  conditional  expectation  with  respect  to  every  a-field  in  the  fil¬ 
tration,  but  for  other  connections  this  may  be  as  difficult  as  solving  a  system  of  nonlinear  partial 
differential  equations,  as  we  shall  now  see. 
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4  Taylor  Approximation  of  a  Gommo-Mortingale  Dirichlet  Problem 
4.1  Condition  for  the  Intrinsic  Location  Parameter 

Consider  a  diffusion  process  {X^,  0  <t  <  5}  on  a  p-dimensionol  manifold  N  with  generator  ^ 
where  A  is  the  Laplace-Beltrami  operator  associated  with  the  diffusion  variance,  and  ^  is  a  vector  field, 
as  in  (7).  The  coordinate-free  construction  of  the  diffusion  X,  given  a  Wiener  process  Won  ,  uses 
the  linear  or  orthonormal  frame  bundle:  see  Elworthy  [1 3]  p.  252.  We  suppose  Xq  =  x^e  N  . 

Also  suppose  (M,  h)  is  a  Riemannian  manifold,  with  Levi-Civita  connection  F,  and  \\t:N  —>M  is  a 
2 

C  map.  The  case  of  particular  interest  is  when  M  =  N,  i|/  =  identity,  and  the  metric  on  N  is  a 
"generalized  inverse"  to  a  •  a  in  the  sense  of  (4).  The  general  case  of  i|/:N  ^  M  is  needed  in  the 
context  of  nonlinear  filtering:  see  Darling  [8]. 

Following  Emery,  and  Mokobodzki  [1  6],  we  assert  the  following: 

W 

Axiom  B:  Any  intrinsic  location  parameter  for  \\f  (X^)  should  be  the  initial  value  of  an  {3^  } - 

2  - 

adapted  hi  F -martingale  {V^,  0  <f  <  5}  on  M,  with  terminal  value  Vg  =  i|/  (Xg)  . 

This  need  not  be  unique,  but  we  will  specify  a  particular  choice  below.  In  the  case  where  a  •  a  (x) 
does  not  depend  on  x,  then  the  local  connector  F,  given  by  (5),  is  zero,  and  Vq  is  simply  E  [\|/  (Xg)  ]  . 
However  our  assertion  is  that,  when  F  is  not  the  Euclidean  connection,  the  right  measure  of  location 
is  Vq,  and  not  E  [\|/  (Xg)  ]  .  We  begin  by  indicating  why  an  exact  determination  of  Vq  is  not  compu¬ 
tationally  feasible  in  general. 


4.2  Relationship  with  Harmonic  Mappings 

For  simplicity  of  exposition,  let  us  assume  that  there  are  diffeomorphisms  (p:N  — ^  R^  and  (p:M 
which  induce  global  coordinate  systems  {x\  ...,x^}  for  N  and  {y\...,y^}  for  M,  respectively.  By 
abuse  of  notation,  we  will  usually  neglect  the  distinction  between  xe  N  and  (p  (x)  e  R^ ,  and  write  x 
for  both.  F  (x)  ( (a  •  a)  (x) )  e  is  given  by  (2)  and  (5),  and  the  local  connector 

F  (y)  e  L  (T^R’^  ®  T^R^;T^R'^)  comes  from  the  Levi-Civita  connection  for  (M,  h)  . 

In  order  to  find  {V^}  ,  we  need  to  construct  an  auxiliary  adapted  process  ,  with  values  in 

R^  j ,  such  that  the  processes  {X^}  and  {  ( V^,  Z^)  }  satisfy  the  following  system  of  forward- 
backwards  SDE: 


X^  =  XQ  +  j*b(X^)ds  +  j‘a(XpdW^,  0<f<5; 


V=^v 


(><8)  -  •  Zp  ds,  0  <  f  <  5  . 


We  also  require  that 


'I 

L  I,  I 


(Vp(Z^-Zp"ds  <-. 


6 


Taylor  Approximation  of  a  Gamma-Martingale  Dirichlet  Problem 


2  - 

[Equation  ( 1 3)  and  condition  ( 1 4)  together  say  that  V  is  on  H  F  -martingale,  in  the  sense  of  ( 1 0) 
and  (1  1 ).]  Such  systems  are  treated  by  Pordoux  and  Peng  [21  ],  but  existence  and  uniqueness  of  solu¬ 
tions  to  ( 1 3)  ore  outside  the  scope  of  their  theory,  because  the  coefficient  F  ( v)  (z  •  z)  is  not  Lipschitz 
in  z. 

2 

However  consider  the  second  fundamental  form  Vd(|)  ofo  C  mapping  Recall  that 

Vd(|)(x)  e  L{T^N  may  be  expressed  in  local  coordinates  by: 

Vd(|)  (x)  (v  ®  w)  =  D^(|)  (x)  (v  ®  w)  -  D(|)  (x)  F  (x)  (v  ®  w)  H-  f  (y)  (D(|)  (x)  v  ®  D(|)  (x)  w)  (15) 

for  (v,  w)  e  ^  '  y  =  ^  ^  Consider  o  system  of  quasilineor  parabolic 

PDE  (a  "heat  equation  with  drift"  for  harmonic  mappings  -  see  Eells  and  Lemaire  [1  1],  [12])  consist¬ 
ing  of  a  suitably  differentiable  family  of  mappings  {u  (f,  .)  :N  M}  ,  for  f  e  [0,  5]  ,  such  that 

^  =  du  ■  ^  +  \vclu  {a  ■  a)  ,  0<t<d ,  (16) 

at  2 

u(0,  .)  =  i|/.  (17) 

For  xe  N ,  the  right  side  of  ( 1 6)  is  du  (f, .)  •  ^  (x)  h-  ^  Vdu  (f,  .)  (a  •  c  (x) )  e  (t  x)  ^  •  f'oHowing 
the  approach  of  Pordoux  and  Peng  [21],  Ito's  formula  shows  that 

Vj  =  u(5-f, Xp  e  M,  (18) 

Z(f)  =  du{d-t,X^)  •a{X^)  e  (19) 

solves  (13).  In  particular  u  (5,  Xq)  =  Vq  .  (A  similar  idea  was  used  by  Thalmoier  [27].) 

4.2. a  Comments  on  the  Local  Solvability  of  (1 6)  -  (1 7) 

Recall  that  the  energy  density  ot  \\t:N  M  is  given  by 

e  (v|/)  (x)  =  ^||di|/  ®  di|/  (a  •  a)  ||^  \  ('l^  M  (x)  (a  •  a  (x) )  .  (20) 

P.y 

Note,  incidentally,  that  this  formula  still  makes  sense  when  a  •  a  is  degenerate.  In  the  case  where 
a  •  a  is  non-degenerate  and  smooth,  ^  =  0,  and  je  (\|/)d(vol^)  <  ,  the  inverse  function  theorem 

method  of  Hamilton  [17],  page  122,  suffices  to  show  existence  of  o  unique  smooth  solution  to  (1 6)  - 
( 1 7)  when  5  >  0  is  sufficiently  small.  For  o  more  detailed  account  of  the  properties  of  the  solutions 
when  dim  (N)  =  2,  see  Struwe  [26],  pages  221  -  235.  Whereas  Eells  and  Sampson  [10]  showed  the 
existence  of  a  unique  global  solution  when  (M,  h)  has  non-positive  curvature,  Chen  and  Ding  [4] 
showed  that  in  certain  other  coses  blow-up  of  solutions  is  inevitable.  The  case  where  a  •  a  is  degener¬ 
ate  appears  not  to  hove  been  studied  in  the  literature  of  variational  calculus,  and  indeed  is  not  within 
the  scope  of  the  classical  PDE  theory  of  Lodyzenskaja,  Solonnikov,  and  Ural'ceva  [20].  A  probabilistic 
construction  of  o  solution,  which  may  or  may  not  generalize  to  the  cose  where  a  •  a  is  degenerate. 
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will  appear  in  Thalmaier  [28].  Work  by  other  authors,  using  Hormander  conditions  on  the  generator 
^  H-  ,  is  in  progress.  For  now  we  shall  merely  assume: 

Hypothesis  I  Assume  conditions  on  a  •  a ,  i|/,  and  h  sufficient  to  ensure  existence  and  unique¬ 
ness  of  a  solution  {u  (t,  .)  —>  M,0  <t  <  ,  for  some  5^  >  0 . 

4.3  Definition:  the  Intrinsic  Location  Parameter 

For  0  <  5  <  ,  the  intrinsic  location  parameter  of  \|/  (Xg)  is  defined  to  be  u  (5,  Xq)  ,  where  Xq  =  Xq  . 

This  depends  upon  the  generator  ^A,  given  in  (7),  where  A  may  be  degenerate;  on  the  mapping 
\|/:N  —>M;  and  on  the  metric  h  forM.  It  is  precisely  the  initial  value  of  on  {3^  }  -adapted  H  F- 
martingole  on  M,  with  terminal  value  Vg  =  i|/  (Xg)  .  However  by  using  the  solution  of  the  PDE,  we 
force  the  intrinsic  location  parameter  to  be  unique,  and  to  have  some  regularity  os  o  function  of  x^  . 

The  difficulty  with  Definition  4.3  is  that,  in  filtering  applications,  it  is  not  feasible  to  compute  solutions 
to  (16)  and  (17)  in  real  time.  Instead  we  compute  on  approximation,  os  we  now  describe. 


4.4  A  Parametrized  Family  of  Heat  Flows 

y 

Consider  o  parametrized  family  {u  ,0<Y<1}  of  equations  of  the  type  (1 6),  namely 


dt 


du^  ■  ^  +  ^Vdu'^  (ct  •  a)  ,  0  <  f  <  5 , 


(21) 


u'^(0,  .)  =  \|/.  (22) 

Note  that  the  cose  y  =  1  gives  the  system  ( 1 6),  while  the  cose  y  =  0  gives  u^  (t,  x)  =  \|/  ((|)j  (x) )  , 
where  { (|)j,  f  >  0}  is  the  flow  of  the  vector  field 

In  o  personal  communication,  Etienne  Pordoux  has  indicated  the  possibility  of  o  probabilistic  con- 
struction,  involving  the  system  of  FBSDE  (46)  and  (55),  of  o  unique  family  of  solutions  {u  }  for  suffi¬ 
ciently  small  y  >  0 ,  and  for  small  time  5  >  0 ,  based  on  the  results  of  Darling  [6]  and  methods  of 
Pardoux  and  Peng  [21].  For  now,  it  will  suffice  to  replace  Hypothesis  I  by  the  following: 

Hypothesis  II  Assume  conditions  on  a  •  a,  \|/,  and  h  sufficient  to  ensure  existence  of  5^  >  0  and 

2  Y 

y j  >  0  such  that  there  is  a  unique  C  mapping  (y,  f,  x)  u  {t,x)  from 
[0,  Yj]  X  [0,  5j]  xN  to  M  satisfying  (2 1 )  and  (22)  for  each  y  e  [0,  y^]  . 

4.4.0  Notation 

For  any  vector  field  ^  on  N,  and  any  differentiable  mop  (():N  P  into  o  manifold  P,  the  "push-for¬ 
ward"  (|)*^  takes  the  value  d(|)  •  ^  (x)  e  T^P  at  y  =  (|)  (x)  e  P;  likewise  (|)*  (^  ®  C,')  =  (|)*^  ®  (|)*^' . 

We  must  also  assume  for  the  following  theorem  that  we  hove  chosen  o  generalized  inverse 
g:T*N  TN  to  a  •  a,  in  the  sense  of  (3),  so  that  we  may  construct  o  canonical  sub-Riemonnion 
connection  V°  for  (.[),  with  respect  tog. 

We  now  state  the  first  result,  which  will  later  be  subsumed  by  Theorem  4.7. 
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4.5  Theorem  (PDE  Version) 

Assume  Hypothesis  II,  and  that  0  <  5  <  5,  .  Then,  in  the  tangent  space  T  ,  ,  M , 

Vdv(xj)  ((|.5.n  )  +I|i.{v"d4.5(x„)  (nj)  -J*((|.5_,).v"d4.,(x„)dn,) 

j -  («, 

‘  y  =  0 

where  { (|)j,  t  >  0}  is  the  flow  of  the  vector  field  (xq)  ,  and 

t  s)  dseT^N®T^N.  (24) 

1  JQ  i  Xj  Xg  Xg 

In  the  special  case  where  M  =  N ,  i|/  =  identity ,  and  h  =  g ,  the  right  side  of  (23)  simplifies  to  the 
part  in  parentheses  {...}. 

4.5.  a  Definition 

The  expression  (23)  is  called  the  approximate  intrinsic  location  parameter  in  the  tangent  space 
denoted  . 

4.5.  b  Remark:  How  the  Formula  is  Useful 

First  we  solve  the  ODE  for  the  flow  {(|)j,  0  <  t  <  8}  of  the  vector  field  compute  3u  (5,  Xq)  /dj  at 
Y  =  0  using  (23)  (or  rather,  using  the  local  coordinate  version  (32)),  then  use  the  exponential  mop  to 
project  the  approximate  location  parameter  on  to  N,  giving 

®’^Pv(x5) 't|o‘^^(S,Xo)  }eM.  (25) 

t  y  =  0 

Computation  of  the  exponential  mop  likewise  involves  solving  an  ODE,  namely  the  geodesic  flow  on 
M.  In  brief,  we  hove  replaced  the  task  of  solving  o  system  of  PDE  by  the  much  lighter  task  of  solving 
two  ode's  and  performing  an  integration. 

4.6  The  Stochastic  Version 

We  now  prepare  an  alternative  version  of  the  Theorem,  in  terms  of  FBSDE,  in  which  we  give  o  local 

coordinate  expression  for  the  right  side  of  (23).  In  this  context  it  is  natural  to  define  a  new  parameter 
2 

e,  SO  that  Y  =  e  in  (21 ).  Instead  of  X  in  (1 2),  we  consider  o  family  of  diffusion  processes 

£  £  2 
{X,e>0}  on  the  time  interval  [0,  5]  ,  where  X  has  generator  ^  h- e  A/2  .  Likewise  V  in  (1 3)  will 

£  2  —  £  £ 
be  replaced  by  o  family  { V  ,  e  >  0}  of  H  F -martingales,  with  Vg  =  i|/  (Xg)  ,  and 

/y  y  — 

Vq  =  u  (5,  Xq)  .  Note,  incidentally,  that  such  parametrized  families  of  F  -martingales  ore  also 
treated  in  recent  work  of  Arnaudon  and  Tholmoier  [2],  [3]. 

4.6.  a  Generalization  to  the  Cose  of  Random  Initial  Value 

Suppose  that,  instead  of  Xq  =  Xq  e  N  as  in  (12),  we  have  Xq  =  exp^  (Uq)  ,  where  Uq  is  a  zero- 
mean  random  variable  in  T  N ,  independent  of  W,  with  covariance  T.„e  T  N  ®T  N ;  the  last 
expression  means  that,  for  any  pair  of  cotangent  vectors  P,  A,  e  T  N , 


Taylor  Approximation  of  a  Gamma-Martingale  Dirichlet  Problem 


9 


E[(p-Uo)  a-U^)]  =  (p®X)  -En. 


Now  set  up  the  family  of  diffusion  processes  {X  ,  e  >  0}  with  initial  values 


^0  =  •  (26 

-  e  ~  W  VV 

Each  F-martingale  {Vj,0<t<5}  is  now  adapted  to  the  larger  filtration  {3f  }  =  {3^  vaiUo)}  . 

In  particular, 

®W(X5)^0 

is  now  a  random  variable  in  depending  on  Uq  . 

4.6. b  Definition 

In  the  case  of  a  random  initial  value  Xq  as  above,  the  approximate  intrinsic  location  parameter  of 
11/  (Xo)  in  the  tangent  space  T  „  x  M ,  denoted  I  y  [\|/  (Xo)  ]  ,  is  defined  to  be 

^  Y  ^0’  ^0  ^ 


-H-p;(x,)Vo. 


We  will  see  in  Section  6.3  below  that  this  definition  makes  sense.  This  is  the  same  as 


|,E[exp;\^^)u''(5,x/"")] 


and  coincides  with  [\|/  (Xg)  ]  ,  given  by  (23),  in  the  case  where  Eq  =  0  . 

4.6. c  Some  Integral  Formulas 

Given  the  flow  {(|)j,  0  <  f  <  5}  of  the  vector  field  the  derivative  flow  is  given  locally  by 

=  (x^)  G  L{T^N;T^N)  ,  C 

where  =  (|)^  (x)  ,  for  0  <  s  <  5  .  In  local  coordinates,  we  compute  as  a  p  x  p  matrix,  given  by 


X*  =  exp{j*D^(x^)du}  . 


Introduce  the  deterministic  functions 


=  K)  (x^f^dsG  T  N®T  N, 


=  +  eT^N®T^N. 


Note,  incidentally,  that  E  could  be  called  the  intrinsic  variance  parameter  of  X  . 
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4.7  Theorem  (Local  Coordinates,  Random  Initial  Value  Version) 

Under  the  conditions  of  Theorem  4.5,  with  random  initial  value  as  in  Section  4. 6. a,  the  approxi¬ 
mate  intrinsic  location  parameter  I  y  [\|/  (X^)  ]  exists  and  is  equal  to  the  right  side  of  (23),  after 

Xg,  Z.g  O 

redefining 

+  (31) 

In  local  coordinates,  21  ^  [\|/  (X^)  ]  is  given  by 

Xq,  o 

jf^[D\{x^)  (sp  -r(xp  (a-a(x^))]dt  +  D\(x5)  (Eg)  -J^^ix^)  (Eq)  +f  (yg)  (Jz/)  (32) 

d  0 

where  >^5  =  t|/  (xg)  ,  J  =  Di|/  (xg) ,  and  E^  is  given  by  (30). 

4.7.0  Remarks 

•  Theorem  4.7  subsumes  Theorem  4.5,  which  corresponds  to  the  cose  Eq  =  0  . 

•  In  the  special  case  where  M  =  N,  \|/  =  identity,  and  h  =  g,  formula  (32)  reduces  to: 

-r(Xf)  (a-a(x^))]df-Tor(xo)  (E^)  +r(xg)  (Eg)}  .  (33) 

•  In  the  filtering  context  [8],  formulas  (32)  and  (33)  are  of  crucial  importance. 

5  Example  of  Computing  on  Intrinsic  Location  Parameter 

The  following  example  shows  that  Theorem  4.7  leads  to  feasible  and  accurate  calculations. 

5.1  Target  Tracking 

In  target  tracking  applications,  it  is  convenient  to  model  target  acceleration  as  an  Ornstein-Uhlenbeck 

process,  with  the  constraint  that  acceleration  must  be  perpendicular  to  velocity.  Thus 

3  3  2 

(v,a)  e  R  xR  must  satisfy  v  ■  a  =  0,  and  the  trajectory  must  lie  within  a  set  on  which  ||v||  is  con- 

2  6 

stant.  Therefore  we  may  identify  the  state  space  N  with  TS  c  R  ,  since  the  v-component  lies  on  a 

2 

sphere,  and  the  o-component  is  perpendicular  to  v,  and  hence  tangent  to  S  . 

Within  a  Cartesian  frame,  X  is  a  process  in  R^  with  components  V  (velocity)  and  A  (acceleration),  and 
the  equations  of  motion  take  the  nonlinear  form: 


df  H-  ®  (34) 

yP{V)dW{t)_ 

Here  the  square  matrix  consists  of  four  3x3  matrices.  A,  and  yare  constants,  W  is  a  three-dimen¬ 
sional  Wiener  process,  and  if  x^  =  (v^,  a^)  , 

p(x)=||o||V||v||", 


dV  _  ®3x3  ^3  V 

dA]  -p(X)/3 -AP(V)  [a 


(35) 
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P(v)  L(R^R^)  ,  (36) 

iivir 

3 

Note  that  P  (v)  is  precisely  the  projection  onto  the  orthogonal  complement  of  v  in  R  ,  and  p  (x)  has 
been  chosen  so  that  d  (V  •  A)  =  0. 

5.2  Geometry  of  the  State  Space 

2 

The  diffusion  variance  metric  (2)  is  degenerate  here;  noting  that  P  =  P,  we  find 


®3x3  ®3x3 


a  =  a  0=  (37) 

03x3YP(-)_ 

—2  6 

The  rescaled  Euclidean  metric  g  =  Y  /e  on  R  is  a  generalized  inverse  to  a  in  the  sense  of  (3),  since 
2  6 

P  =  P  .  We  break  down  a  tangent  vector  ^  to  R  into  two  3-dimensional  components  and  . 

2 

The  constancy  of  ||v||  implies  that 

DP(v)\  =  +  •  (38) 


Referring  to  formula  (5)  for  the  local  connector  F  (x)  , 


D<g(C)lg(?))(il)  =  +  (C.  ?.il)  e  T^NxT^NxT^N. 

Y  llvll 


Taking  first  and  second  derivatives  of  the  constraint  v  •  o  =  0 ,  we  find  that 

^^vH-C^o  =  0,  +il^C  =  0- 

Tq  '3y  /  Iq  Ty  ly  Tq 

Using  the  last  identity,  we  obtain  from  (5)  the  formula 


r(x)  (C®q)  =  S(C®q)  =  + 


_r  r"'’ 


In  orderto  compute  (32),  note  that,  in  particular. 


r(x)  (a  a(x))  =  ^  V  =  ®  . 


5.3  Derivatives  of  the  Dynamical  System 

It  follows  from  (7),  (34),  and  (41 )  that  the  formula  for  the  intrinsic  vector  field  ^  is: 
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p  (x)  V  -  A,P  (v)  a 


Differentiate  under  the  assumptions  ||v||  is  constant  and  v  ■  a  =  0,  to  obtain 


D^(x)  (Q  =  r  °  '  1  ,  Q^—^- 

[XQ-p-XP-2qJ  llvll' 


Differentiating  (43)  and  using  the  identities  (39), 


D  ^(x)  (ii®Q 


IMI^  + 


5. 3. a  Constraints  in  the  Tangent  Space  of  the  State  Manifold 

Let  us  write  a  symmetric  tensor  xe  T  N  ®T  N  in  3x3  blocks  as  the  matrix 


Xvv  ^vo 
^ov  ^ao 

where  xly  =  •  Replacing  by  Xya-  and  by  Trfx^,^)  ,  etc.,  in  (44),  we  find  that 


d\{x)  ix)  =  ^  ^  ° 

IMI^  Tr(X„„)v+  {Xya  +  Xay^° 


5.4  Ingredients  of  the  Intrinsic  Location  Parameter  Formula 

6  2  6 
Le\  \\f  :N R  be  the  inclusion  of  the  state  space  N  s  TS  into  Euclidean  R  .  Thus  in  formula  (32), 

-  2 

the  local  connector  r(.)  is  zero  on  the  target  manifold,  J  is  the  identity,  and  D  \|/  is  zero.  When  Eq 

is  token  to  be  zero,  the  formula  for  the  approximate  intrinsic  location  parameter  mg  =  [Xg]  for  Xg 
becomes: 


^  ''’r(X^„(f))Vj+  (XyaW  +Xav(^))°t 


where  and  o^  are  the  velocity  and  acceleration  components  of  x^  =  (|)j  (x) ,  for  0  <  f  <  5 ,  and 
and  Xt  given  by  (28)  and  (29).  A  straightforward  integration  scheme  for  calculating  (t[^,  x^, 
at  the  same  time,  using  a  discretization  of  [0,  5]  ,  is: 

i[  =  exp  [D^  (x^)  +  (xj)  ]  }  , 


^  (a  •  a)  (xp  +  Xu  + 
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m. 


m 


+  ■ 
I  2 


Since  the  local  connector  is  zero  on  the  target  manifold,  geodesics  ore  simply  straight  lines,  and 
Xg  +  [Xg]  is  0  suitable  estimate  of  the  mean  position  of  Xg . 

5.5  Simulation  Results 

FIGURE  1  SIMULATIONS  OF  THE  MEAN  OF  AN  SDE,  VERSUS  ITS  APPROXIMATE  ILP 


x-acceleration:  ODE  mean  SDE  +/- 1  s.e. ILP  "o" 


4  3 

We  created  10  simulations  of  the  process  (34),  with  A,  =  0.5  and  y  =  5.2  x  10  ,  on  the  time  inter¬ 
val  [0,  1]  ,  which  was  discretized  into  25  subintervals  for  integration  purposes.  In  each  cose  V  and  A 

'y 

were  initialized  randomly,  with  magnitudes  of  200  m/s  and  50  m/s  ,  respectively.  The  plot  shows  the 
x-component  of  acceleration  (the  other  two  ore  similar):  the  "  +  "  signs  represent  the  mean  of  the  pro- 

4 

cess  (34)  over  10  simulations  ,  with  bonds  showing  plus  and  minus  one  standard  error,  the  solid  line 
is  the  solution  {x^,  0  <f  <  1}  of  the  discretized  ODE  dx^/dt  =  ^  (x^)  ,  and  the  circles  denote  the 
approximate  intrinsic  location  parameter  (ILP).  The  reader  will  note  that  the  ILP  tracks  the  mean  of 
the  process  better  than  the  ODE  does. 
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6  Proofs  of  Theorem  4.7 

The  strategy  of  the  proof  will  be  to  establish  the  formula  (32)  using  Ito  calculus,  and  then  to  show  it  is 
equivalent  to  (23)  using  differential  geometric  methods.  While  this  may  seem  roundabout,  the  impor¬ 
tant  formula  for  applications  is  really  (32);  converting  it  into  (23)  serves  mainly  as  o  check  that  for¬ 
mula  (32)  is  indeed  intrinsic.  It  will  moke  no  difference  if  we  work  in  global  coordinates,  and  identify  N 
with  and  M  with  . 

6.1  Step  I:  Differentiation  of  the  State  Process  with  Respect  to  o  Parameter 

We  consider  o  family  of  diffusion  processes  {X  ,  e  >  0}  on  the  time  interval  [0,  5]  ,  with  initial  values 
£ 

Xq  =  exp^  {eUq)  ;  here  is  o  zero-mean  random  variable  in  N ,  independent  of  W,  with  covori- 

0  g  2  ^ 

anceIl„eT  N®T  N,  and  X  has  generator  ^  h- e  A/2  . 

U  Xg  Xg 

Note  that,  in  local  coordinates,  the  SDE  for  X^  is  not  "dX^  =  b  (X^)  ds  +  ea  (X^)  dW^  ",  because  the 
limiting  case  when  e  =  0  would  then  be  the  ODE  based  on  the  vector  field  ^b'O. ,  which  is  not  the 
some  as  which  is  given  by  (7).  Instead  the  SDE  is 

X"  =  exp^^  (eUo)  +  (X")  ds  -  C  (X")  ds  +  (X")  dW^  (46) 

where  we  use  the  notation 

C(x)  =^r(x)  (a(x)  afx))  .  (47) 

In  the  cose  e  =  0 ,  the  solution  is  deterministic,  namely  {x^,  0  <  f  <  5}  .  Note  that,  in  local  coordi¬ 
nates, 

2 

expx^(eDo)  =  XQH-eUQ-yr(xQ)  (Uq®  Uq)  H-o(e^)  . 

It  is  well  known  that,  if  the  vector  field  ^  and  the  semi-definite  metric  (.|.)  ore  sufficiently  differentio- 

E  2  e  2 

ble,  then  the  stochastic  processes  3X  /3e  and  3  X  /3e  exist  and  satisfy  the  following  SDEs: 


?  =  +  Is  +  f  a(X;)dW, 


-2e/  CCXjjds  +  ef  Do(X;)  dW,  +  0(e")  ; 


2  g  /  £  £\  f  ^  £ 

O  X,  rt  0  P  rt  p  d  X  fi  P 

— *  =  -r(x)  +j^D^^(X:)^^^®g^^Jds  +  j^D^(X:)  — Hds-2f^C(X:)ds  (49) 
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where  0(e)  denotes  terms  of  order  e.  Define 


Now  (48)  and  (49)  give: 


e  =  o'  *  ^ae^ 


dAj  =  AjAjdf  +  a (Xj) dWj,  Aq  =  U^; 


dA^ ^  =  [d\  (Xj)  (  Aj  ®  Aj)  +  AjAf  ’  -  2C  (xj)  ]  df  +  2Da  (x^)  ( A^)  d W, ,  (52) 

A^^  =  -r(xo)  (Do®Uo)  , 

where  Aj  =  D^(Xj)  .Let  {T*,0<s,t<5}  be  the  two-porometer  semigroup  of  deterministic  matrices 
given  by  (28),  so  that 


"  b  s  ,  r  r  t 

—  =  — XA'T  =  TT 

3^  t  t'  s  r  s 

Then  (51)  becomes  d(x|*Aj)  =  i^a  (x^)  dW^,  which  has  a  Gaussian  solution 

^  =  ^*0^0  +  ^*0 j*  =  ^1^0  +  Jo •  (53) 

Likewise  (52)  gives  d(T[*Aj^^)  =  [D^^(Xj)  {A^®  A^)  -  2^  (x^)  ]  dt  +  2x|*Da  (x^)  (A^)  dW^,  whose 

solution  is 

A^^  =  -^!)r(V  (A,®A^)  -2C(x3)]ds  +  2Da(x^)  (A^dWJ  .(54) 

6.2  Step  II:  Differentiation  of  the  Gommo-Martingole  with  Respect  to  a  Parameter 

£  £  £ 
Consider  the  pair  of  processes  (V  ,Z  )  obtained  from  (18)  and  (19),  where  u  is  replaced  by  u  .As 

£  £ 

in  the  case  where  e  =  1 ,  ( V  ,  Z  )  gives  an  adapted  solution  to  the  backwards  equation  correspond¬ 
ing  to  (13),  namely 

vj  =  v(X5)-)^z;dw,+i(f(v;)  (z^zj)*. 

£  £ 

However  the  version  of  ( 1 9)  which  applies  here  is  Z^  =  Du  ( 5  -  f,  X^)  ea  (X^)  ,  so  we  may  replace 
£  £ 

Z^  by  eZ^ ,  and  the  equation  becomes 
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V*  =  ¥  (Xg)  -  ej'z^^cfW^  +  J^f  ( ds .  (55) 

£  £  £  2,  £  2,  2.  £  2, 

By  the  regularity  of  u  ,  it  follows  that  dV  /de,  dZ  /de,  d  V  /de  ,  and  d  Z  /de  exist,  and  satisfy 

the  following  equations: 


^  =  Dw  (Xg)  I  ^ ^  -  j;z:dW^  -  ej^  ^^dW^  +  e  1/  ( Vp  (Z  •  Z^  ds+  O  (e  )  ;  (56) 


e'3Xg  2  gfax^g  axn  gdz"  e  e  e 

>=Dv(X5)  +Dv(X5).  -2  jj‘dW,+  ^r(V.)(Z,.Z,)ds+0(e).,S7, 

0  ~  w  vv 

Note  also  that  Vj  =yg  =  i|/(xg)  for  all  fe  [0,  5]  .  Take  3f  =5^  v  a  (Ug)  .  By  combining  (53)  and 
(56),  we  see  that,  if  J  =  Di|/  (xg)  , 


Define 


oV.  r  ~VV“|  5  i-t  8 

®,-35  =  eK|3,  ]  z;a(.^)dW^: 

e  =  0 


7O  ;  8  ,  , 

Z  =  Jx  a  {x  )  . 

S  S 


0(2) 

*  32 

de 


r  (2)  1  rS-  00 

E  '  +  D  V  (Xg)  (Ag  ®  Ag)  +  f  r  (yg)  (Z^  •  Z^  ds  3f 


From  (54)  and  (61 )  we  obtain: 


0*^’  =  -jTgr(xg)(Ug®Uo)+JEr{j%^^[D^^(x^)(A^®A^)-2C(x^)]ds}  Ug 


+  E  {dV(x8)  (Ag®Ag)  +  jgr(yg)  (z“-z°)ds}  Ug  . 


The  expected  value  of  o  quadratic  form  q  Aq  in  on  (p,  L)  random  vector  q  is  easily  computed  to 
be  +  p^Ap  .  In  this  case. 


q  =  A^~Np(TgUg,X3), 


where  %  is  given  by  (29),  so  we  obtain 


32)  _  ,A 


-JTgr(xg)  (Do®Ug)  +  J  [ D (x^  (^3)  +  (x {x^U ^  ®  x^U -2C(x^)]ds 
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+  d\  (xg)  (Xg)  +  d\  (xg)  (TqUo  ®  XqUo)  +  r  (yg)  UX^/)  ■  (6i 

6.3  Step  III:  A  Taylor  Expansion  Using  the  Exponential  Mop 

Let  yg  =  i|/ (xg)  ,  and  define  P(e)  ^Ej^exp^V^J  .  Referring  to  (27),  we  ore  seeking  3p/5e^|g  _Q.  It 
follows  immediately  from  the  geodesic  equation  that 

D(expy)  \y)  (w)  =  w ,  (exp^)  ^(y)(v®w)  =f(y)(v®w)  (63 

It  follows  from  (58)  and  (62)  that 


./E  ,  8, ,  e  ,-.(2)  ,2. 

Vq  =  yg  +  eJ'^gL/o  +  -©0  +  o  (e  )  . 


A  Taylor  expansion  based  on  (63)  gives 


-i,,e  .  1 


e>^Py,^o  =  ^0-y8+5^(y8)  ((^0-ys)  ®  (^0-^5))  +ol  ^0-^8 


=  eJzlu^  +  I  {0*^’  +  r  (yg)  (JtJUq  ®  JToUq)  }  +  o  (e^)  . 

Taking  expectations,  and  recalling  that  has  mean  zero  and  covariance  ^  N  ®T  N ,  we 

U  U  Xp  Xg 

obtain 


It  follows  that  ^ 


p  (e)  =  ^  {  E  [0^^^  ]  +  r  (yg)  (jxolo  (JTq)  ^J}  +  o  (e^) 


0,  and  hence  that 


=  ^  [d\  (xp  (xp  -  r  (xp  (a  •  a  (xp  )]ds  +  d\  (xg)  (Xg)  +  r  (yg)  Ux/) 

+  jf^%lD\{x^)  (ipds-JTor(xo)  (Ig)  +DV(xg)  (Ig)  +f(yg)  (JEg/)  } 

s  s  2" 

where  =  '^0^0  ^"^0^  •  write 


S  =x  +  , 

s  s ' 


then  the  formula  becomes 
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i  [d\  (xp  (sp  ]  ds  +  D  V  (Xg)  (Sg)  +  r  (yg)  (JSg/)  } 


- \j  {  (x^)  (a  •  a  (xp )  ds  +  T^r  (xq)  (Iq)  }  . 

This  establishes  the  formula  (32).  0 

6.4  Step  IV:  Intrinsic  Version  of  the  Formula 

It  remains  to  prove  that  (23),  with  os  in  (3 1 ),  is  the  intrinsic  version  of  (32).  We  abbreviate  here  by 
writing  i|/  •  (|)^  os  i|/j .  By  definition  of  the  flow  of 

^\|/j(xo)  =  di|/«^((^j(xQ))  , 

and  so,  differentiating  with  respect  to  x,  and  exchanging  the  order  of  differentiation. 


=  D(|Vf(Xo) 


or,  by  analogy  with  (28),  taking  0^  =  ^^t-s  (x^)  •  D  (^l^t  • 


:^)  (X  )  =  Dv(Xj)t' 


dx  f  ae  f 

(xp  x^ ;  =  Dv  (xp  •  (x^)  x^  . 

Since  x^Xq  =  Xq,  we  have  G^Xq  =  Di|/ (xg)x^XQ  =  0q,  which  upon  differentiation  yields 


h‘ >0  = -0.  5  = 


which  gives 

^  =-T,D^(x,);^^  =-0,D^(xp.  ( 

A  further  differentiation  of  (65)  when  i|/  is  the  identity  yields 

(x^)  ]  (v  ®  w)  =  d\  (xp  (XqV  ®  XqW)  +  (xp  (xq)  (v  ®  w)  .  ( 

Combining  (67)  and  (68),  we  have 

^(efD\(xo))  =  -efD^(xpD\(xo)  +ef{D^^(xp  (x*  (.)  ®Xo(.))  +D^(xpD\(xo)} 


|^(efD\(xo))  (V®w)  =  0,V^(xp  (XoV®Xow)  . 


The  formula  (15)  for  V°d(|)j  con  be  written  as 
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V°d(^^(Xo)  (v®w)  =  d\(Xo)  (v®w)  -Tor(xo)  (v®w)  +r(x^)  (V®T*w) 


It  is  clear  that 

|(efx;,r(xo))  =^(eSr(xo))  =  o. 

Hence  from  (69)  and  (70)  it  follows  that 

|r{(V8_f)*V°d(^j(xo)}  (v®w)  =  |-{0f[D\(xo)  (v®w)  +T{x^)  (v®ToW)]} 


=  (xp  (Tqv  ®  Tqw)  +  i { efr  (x^)  ( ^  ®  x^w) } . 


The  last  term  in  (71)  can  be  written,  using  (66),  as 

^{efr(xj)}  (TqV®i:qW)  +0jr(Xj)(^D^(xpTQV®TQW  +  i:QV®D^(xpT*wj, 

for  V  ®  w  e  T  N  ®T  N  .  We  will  replace  v  ®  w  by 
^0  ^0 

n* -  ^  =  ^0  +  lo^  ^  ®  ^x,^ 

t  t  T 

where  Z^  =  x^  +  XqEq  (Tq)  .  Observe  that 

^ont(4)^  =  S,eT  N®T  N. 

Moreover  from  (74)  and  (66),  it  is  easily  checked  that 


-'=  (c-a)  (xp +D^(x^)S^  +  S^{D^(x^)}  . 


It  follows  from  (71)-  (75)  that 

(  dn  ^ 

|{(V8_r)*V°dc^,(x,)  (n^)}  -  (¥8_r)*V°dc^r(xo)  U  M  =  |{  (^8-*)  *V°d(^,(x,) }  (H^) 


e^D\{x^)  +|-{0fr(xp}l(Ep  +0fr(Xj)  [D^(xpE^+E^(D^(xp)^] 


(xp  +  ^{efr (xp }  (Sj)  +  0fr (x^) I  (a  •  a)  (xp 


0fD^^(xp  (sp  -0fr(x^)  ((a-a)  (x^)  +^{0fr(xp  (S^} 
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Since  V°d(|)Q  =  0,  it  follows  upon  integration  from  0  to  5  that  in  , 

V*V°d(^g(xo)  (Hg)  -j®(Vg_P*(V°d(^^(xo))dn^  = 

Dv  (Xg)  { [d\  (xp  (Sj)  -  r  (x^)  {{a  a)  (x^) )  ]  dt  +  r  (xg)  (Eg)  -  x^r  (xq)  (Eq)  }  .  (76) 

However  the  formula  ( 1 5)  for  Vdi|/  (xg)  (v  ®  w)  can  be  written  as 

D^\|/  (xg)  (v  ®  w)  -  Jr  (xg)  (v  ®  w)  +  r  (yg)  (Jv  ®  Jw)  .  (77) 

where  )^g  =  (xg)  ,  and  J  =  D\\i  (xg)  .  We  take  v  ®  w  =  Eg  e  N  ®T^  N ,  and  add  (76)  and  (77): 

J  [®xf  [d\  (xp  (E^  -  r  (xp  {{a  a)  (x^ )  ]  df  +  D  V  (xg)  (Sg)  -  Jxjr  (x^)  (E^)  +  f  (yg)  (Jz/) 

=  vdv(xg)  (((^g)^ng) +v*{v°d(^g(xo)  (Hg) -j®((^g  ^)  jv°d(^^(xo))dn^}  . 

The  equivalence  of  (23)  and  (32)  is  established,  completing  the  proofs  of  Theorems  4.5  and  4.7.  0 

7  The  Canonical  Sub-Riemannian  Connection 

The  purpose  of  this  section  is  to  present  a  global  geometric  construction  of  a  torsion-free  connection 

2 

V°  on  the  tangent  bundle  TN  which  preserves,  in  some  sense,  a  C  semi-definite  metric  (.|.)  on  the 
cotangent  bundle  T*N  induced  by  a  section  a  of  Horn  (R^;TN)  of  constant  rank.  In  other  words, 
we  assume  that  there  exists  a  rank  r  vector  bundle  E  —>  N ,  a  sub-bundle  of  the  tangent  bundle,  such 
that  =  range  (a  (x) )  c  T^N  for  all  x  g  N  . 

Given  such  a  section  c,  we  obtain  a  vector  bundle  morphism  a:T*N  TN  by  the  formula 

a  (x)  =  a  (x)  •  i  •  a  (x)  * ,  X  G  N ,  (78) 

where  i:  (R  )  R  is  the  canonical  isomorphism  induced  by  the  Euclidean  inner  product.  The  rela¬ 
tion  between  a  and  (.|.)  is  that,  omitting  x, 

p  •  a  ( A.)  =  (plA) ,  Vjj.  G  T*N  .  (79) 


7.1  Lemma 

Under  the  constant-rank  assumption,  any  Riemannian  metric  on  N  induces  an  orthogonal  splitting  of 
the  cotangent  bundle  of  the  form 


T^*N  =  Ker(a(x))  ©F^,  (80) 

where  E  ^  N  is  a  rank  r  sub-bundle  of  the  cotangent  bundle  on  which  (.|.)  is  non-degenerate.  There 
exists  a  vector  bundle  isomorphism  a°:T*N  TN  such  that  E^  =  Ker  (a°  (x)  -  a  (x) ) ,  and 
a°  (x)  ^  is  a  generalized  inverse  to  a  (x) ,  in  the  sense  that 
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a  (x)  •  a°  (x)  ^  •  a  (x)  =  a  (x)  . 

Proof:  For  any  matrix  A,  range  (A)  =  range  (AA^)  ,  and  so  range  (a  (x))  =  range  (a  (x))  .  It  fol¬ 
lows  that 

dim  Ker(a(x))  =  p-dimE^  =  p-r. 

Given  a  Riemannian  metric  g  on  N  (which  always  exists),  let  (.|.)°  be  the  dual  metric  on  the  cotan¬ 
gent  bundle.  Define 

F^=  {0G  T/N  :  =  0  VA,e  Ker(a(x))}  . 

We  omit  the  proof  that  F  ^  N  is  a  vector  bundle.  Since  dim  Kerf  a  (x) )  =  p -r,  it  follows  that 
dimF^  =  r;  since  the  rank  of  (.|.)  is  r,  we  see  that  (6|0)  >  0  for  all  non-zero  0  e  F^ .  This  shows  that 
(.|.)  is  non-degenerate  on  the  sub-bundle  F  N  . 

Now  (80)  results  from  the  orthogonal  decomposition  of  with  respect  to  (.|.)°  .  Hence  an  arbi¬ 

trary  X  e  T^*N  can  be  decomposed  as  A,  =  Aq  ©  Aj ,  with  Aq  e  Ker  (a  (x) )  and  Aj  e  F^  .  The  metric 
(.|.)°  induces  a  vector  bundle  isomorphism  ^:T*N  TN ,  namely 

^•p(A)  =<p|A)°,  Vpe  T*N. 

Now  define  a°:T*N  —>TN  by 

a°(A)  =P(Ao)  H-a(Aj)  . 

It  is  clearly  linear,  and  a  vector  bundle  morphism.  Since  P  is  injective, 

{Ae  T/N  :  a°(A)  =a(A)}  =  {A^  ©  A^  e  T/N  :  p  (A^)  =  0}  =  F^, 

which  shows  that  F^  =  Ker  (a°  (x)  -  a  (x) )  .To  show  a°  is  an  isomorphism,  it  suffices  to  show  that 
a°  ( A)  5^  0  whenever  A  5^  0  .  When  Aq  0 ,  non-degeneracy  of  (.|.)°  implies  that 

AQ-a°(A)  =  <Aq|Aq)°  +  (AqIAj)  =  <Aq|Aq)°^0. 

On  the  other  hand  when  Aq  =  0 ,  and  Aj  5^  0 ,  non-degeneracy  of  (.|.)  on  F  N  implies  that 

Aj-a°(A)  =  A^-a(A^)  =  <ApA^)?::0. 

Hence  a°:T*N  TN  is  a  vector  bundle  isomorphism  as  claimed.  The  generalized  inverse  property 
follows  from  the  fact  that  a  (x)  •  a°  (x)  ^‘afxffA)  =  A^ .  0 

7.2  Proposition 

Suppose  o  is  a  constant-rank  section  of  Horn  (R^;TN) ,  inducing  a  semi-definite  metric  (.|.)  on  T*N 
and  a  vector  bundle  morphism  a:T*N  TN  as  in  (78)  and  (79).  Suppose  furthermore  that 
a°:T*N  ^  TN  is  a  vector  bundle  isomorphism  such  that  a  (x)  •  a°  (x)  ^  •  a  (x)  =  a  (x) ,  as  in 
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Lemma  7.  /.  Then  TN  admits  a  canonical  sub-Riemannian  connection  V°  for  (.|.),  with  respect  to 
a° ,  which  is  torsion-free,  and  such  that  the  dual  connection  V  preserves  (.|.)  in  the  following  sense: 
for  vector  fields  V  in  the  range  of  a,  and  for  I -forms  6,  X  which  lie  in  the  sub-bundle  F  =  Ker  ( a°  -  a) , 

V<0|X)  =  (VvGlA,)  +  (eiVyX) .  (81) 

[Here  V^d  •  W  =  Z  (0  •  W)  -  6  •  V°^W.J  For  any  l-forms  0,  |a,  X,  and  corresponding  vector  fields 


Y=a°{e)  ,  Z  =  a°{[i)  ,  W=a°{X)  , 


the  formula  for  V°  is: 


^  =  ^  {  Vai^)  +  Wim  -  Z{Q\X)  +  X  •  [Z,  Y]  +  ^  •  [Y,  W]  -  0  •  [W,  Z]  }  .  (82) 

7.2. a  Expression  in  Local  Coordinates 

Take  local  coordinates  for  N,  so  that  a°  (x)  ^  is  represented  by  a  matrix  (g/^)  ,  and  a  (x)  by  a 
ik 

matrix  (a  )  ,  where  by  Lemma  1  FI, 

Ejk  rm  jm 

a  gi^^a  =  a  . 


Take  Y=d/dx.,  W  =  d/dx. ,  and  Z  =  d/dx^  in  (82),  so  that  p,  =  becomes 

L„s  1 r  3  ,  rs  3  rs  3  rs 

r,7  9si<=  22-^37 (9/r“  Ssfc) +37(9/r«  9sfc)-^(9,>a  9,y)  >  • 

r,  s  i 


When  (.|.)  is  non-degenerate,  then  529/r^^^  “  reduces  to  the  standard  formula  for  the 

Levi-Civita  connection  forg,  namely  ^ 


lldgjkjQik  ^9ii 


^  2  3x.  3x.  3xJ  ■ 

S  ^  !  j 


7.2.  b  Remark 

A  similar  construction  appears  in  formula  (2.2)  of  Strichartz  [25],  where  he  cites  unpublished  work  of 
N.  C.  Gunther. 

7.2. C  Proof  of  Proposition  7.2 

First  we  check  that  the  formula  (82)  defines  a  connection.  The  R-bilinearity  of  ( Y,  W)  V°yW  is 
immediate.  To  prove  that  Y°yfY/  =  fV°YW  +  (Yf)  W  for  all  f  e  C  (N)  ,  we  replace  W  by  fW  and  X 
by  fX  on  the  right  side  of  (82),  and  the  required  identity  holds.  Verification  that  (82)  is  torsion-free  is 
likewise  a  straightforward  calculation. 

Next  we  shall  verify  (81 ).  By  (79),  and  the  definition  of  duality,  taking  V  =  a  (0)  , 


)t-V°ya(A,)  =  V(p  •  a(A,))  -  Vyp  •  a(A,)  =  V<A,|)t)-<A,|Vv)t) . 


Switch  )j.  and  X  in  the  last  expression  to  obtain: 
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A,-V°^a()t)  =  V<A,|)t)-<|t|VvA,). 

Taking  U  =  a  ()t)  and  S  =  a  (A.)  ,  we  see  that 

V<A|^)  -  <A|Vv^)  -  (^IVyA)  =  ^  •  V°yS  +  A  •  V\U  -  V<A|^) .  (84) 

In  terms  of  the  splitting  T^*N  =  Ker  (a  (x) )  ©  of  Lemma  7. 1 ,  we  may  write  9  =  6q  +  6^ ,  etc.,  and 
we  find  that,  if  oc  (0)  ,  then  V  =  a  (0^)  =  a°  (0^)  ,  etc.  It  follows  from  (82)  that 

fii- V°yS  =  ^{V<AilPi)  +  S<^iiei)-U<ei|Ai)  +  Ai-  [U,V]  +^ij-  [V,S]  -e^-  [S,U]}; 

Ai  •  ^  {  V<AilPi)  +  U<Ai|0i)  -  S<eil^i)  +  Pi  •  [S,  V]  +  Aj  •  [V,  U]  -  01  •  [U,  S]  }  . 

To  prove  (81),  we  can  assume  Aq  =  Pg  “  now  the  right  side  of  (84)  becomes 

Pi-V°^S  +  Ai- V°^U-V<AilPi)  =  0, 


os  desired.  0 

8  Future  Directions 

We  would  like  to  find  out  under  what  conditions  on  a  •  a,  \|/,  and  h  the  system  of  PDE  (1 6)  -  (1 7) 
has  a  unique  solution  for  small  time,  other  than  the  well-known  cose  where  a  •  a  is  non-degenerate, 
and  ^  =  0;  likewise  for  the  parametrized  family  (21 )- (22).  It  is  likely  that  the  conditions  will  involve 
the  energy  of  the  composite  maps  {\|/ •  (|)i,  0  <  t  <  5}  .  Both  stochastic  and  geometric  methods  should 
be  considered.  Another  valuable  project  would  be  to  derive  bounds  on  the  error  of  approximation 
involved  in  the  linearization  used  in  Theorem  4.7.  This  is  likely  to  involve  the  curvature  under  the  diffu¬ 
sion  variance  semi-definite  metric  -  see  Darling  [9]. 
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